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The objective of this research is to develop turbulence models to predict the flow and heat 
transfer fields dominated by the curvature effect such as those encountered in turbine cascades and 
turn-around ducts. 

A Navier-Stokes code has been developed using an explicit Runge-Kutta method with a two 
layer k-e/ARSM (Algebraic Reynolds Stress Model), Chien’s Low Reynolds Number (LRN) k-e 
model and Coakley’s LRN q- « model. The near wall pressure strain correlation term was included 
in the ARSM. The formulation is applied to Favre- averaged N-S equations and no thin-layer 
approximations are made in either the mean flow or turbulence transport equations. Anisotropic 
scaling of artificial dissipation terms was used. Locally variable timestep was also used to improve 
convergence. Detailed comparisons were made between computations and data measured in a turbine 
cascade by Arts et al. at Von Karman Institute. The surface pressure distributions and wake profiles 
were predicted well by all the models. The blade heat transfer is predicted well by k-e/ARSM 
model, as well as- the k- e model. It’s found that the onset of boundary layer transition on both 
surfaces is highly dependent upon the level of local freestream turbulence intensity, which is strongly 
influenced by the streamline curvature. 

Detailed computation of the flow in the turn around duct has been carried out and validated 
against the data by Monson as well as Sandborn. The computed results at various streamwise locations 
both on the concave and convex sides are compared with flow and turbulence data including the 
separation zone on the inner well. The k-e/ARSM model yielded relatively better results that the 
two-equation turbulence models. A detailed assessment of the turbulence models has been made with 
regard to their applicability to curved flows. 
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OBJECTIVE : 


TO DEVELOP TURBULENCE MODELS TO PREDICT 
FLOW AND HEAT TRANSFER FIELDS IN 
TURBOMACHINERY INCLUDING CURVATURE, 
ROTATION AND HIGH TEMPERATURE EFFECTS 


OUTLINE : 

• INTRODUCTION 

• NUMERICAL TECHNIQUE 

• TURBULENCE MODELS 

• FLOW AND HEAT TRANSFER FIELD IN A HIGH 
MACH NUMBER TRANSONIC TURBINE CASCADE 

• FLOW FIELD IN A TURN-AROUND DUCT 

• CONCLUSIONS 
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GOVERNING EQUATIONS 
(Cartesian) 
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Relative velocities, constant rotation rate about x axis, co. Averaged quantities. 
Energy , tQmt + 2 ' 2 .rothalpy constant along streamlines for inviscid steady state. 



TECHNIQUES 


1 . RK2D code : 

* 2-D Navier-Stokes code. 

Conservative, compressible formulation 

* Favre-Averaged Mean and Turbulence 
equations 

* 4-stage explicit Runge-Kutta scheme 

* 2nd and 4th order artificial dissipation 
(with eigenvalue and local velocity scaling) 

* Coupled with compressible Low- 

Reynolds number K-e model , q-co model , 
ARSM , NLSM ( Nonlinear -stress model) , 
AHFM(Algebraic Heat Flux model) 

* Characteristic boundary conditions, H 
grids (generated by a combined algebraic 
and elliptic method to keep smoothness and 
orthogonality near the wall) 


2. TEXSTAN code 


* 2-D boundary layer code developed by 
Crawford 

* Extension of STAN 5, Patankar- 
Spalding numerical scheme 

* Include 7 differential two-equation 
turbulence models (Jones-Launder, 
Chien,Lam- Bremhorst , etc.) and mixing 
length model 
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RSM ( Reynolds Stress Model) 
(Gibson & Launder 1978) 

Reynolds stresss transport equation : 


duu 


U k -5^ = + D Kj +u j,i } 


3x, 


P U j~ PU; 5 

U : U.U,, + + — -*5 V — =r 


3u : U; 


■i“j“k 


ik 


p jk v a x . 


3uj du. 
-2v — * i 


3x k 3x k 


1,6, , C-. — D-. — P.. + ~ 

’ ij ij u y ij ij 


where e..=|e (Dissipation) 

J 

‘t’ij =< t ) iji +< l ) ij 2 +< l ) iji,w +< t , ij 2 ,w (Pressure-strain 
correlation) 


4 >jji =- c i|(E^u]-^ Wy) (Retum-to-isotropy 
part) 

<t>i j2 = -c 2 (Py -fp^y) (Rapid part) 
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fyjw.l _C l^ u k UmI, k n m^ij~f u i u k n k n j - § u j u k n k n i^n 

(Near-wall term) 

i 3 3 

^ijw,2 = c 2^ < t 1 km,2 n k n "'^ij ~2 < l , ik,2 n k n j ~ 2^jk,2 n k n i^ n 

(Near-wall term) 

f n = k 3/2 /(2.55x n e) (x„ is the distance normal 
to the wall) 

Constants : c x =1.8, c 2 =0.6, cj =0.5, c 2 
=0.3 
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ARSM (Algebraic Reynolds Stress Model 



ARSM assumption: 


c ij“ D ij = ^ (C k _D k )= -ir (P k- e) 


U : U. 


:> c ij- D ij=nr( c k- D k)=4r i(P «- e) 


U;U ; U;Uj 

:> c y _D u = "Tr 1 ^ _D k ) = 4r (P k - £) 


where 


P, =— u u U. • 

k i Ji,J 


P U = - u i u k U j > k _i W J i,k 




and c|) ij2iW as in RSM 
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NLSM (Nonlinear Stress Model") 
(Shih , Zhu & Lumley 1992) 


Reynolds stress : 


+- 


+- 


+- 


u]u] = |k5 ij -v 1 (U. j + U. i ) 

(U k,i U k,j 
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where 


7C = U. -U-- 

»>J J>1 

7 t = U . - U . • 

i,J i.J 
k 2 

v t =C^ 


C„ = 


2/3 


** Aj+rj+cd; 


t = -n 

^ e 


Q = (2Q* j a* j )2 

a ij =(u i, j - u j,i )/2 

^ k c 
r\ = —b 
e 


S = ( 2 S ij S ij )2 
S ij =( U, j + U j , i ) / 2 


Constants: 


C Tl 

C ,2 

C t2 

A x 

a 

A 2 

-4 

13 

-2 

1.25 

0.9 

1000 


1782 












turbulence intensity 



Fig. Validation of the ARSM model: 
Turbulence intensity profile in the flat-plate 
turbulent boundary layer : experiment by 

Klebanoff; computation by ARSM 
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1 80 -riegree TURN AROUND DUCT (TAD 
Geometry & Grid : 



EARLIER RESEARCH ON TAD ROW 


• Measurements: 

*Sandbom (1988), Sandbom and Shin 
(1989) (Water flow, Re=7xio 4 - s*io 5 (Re 
based on duct height and bulk velocity) 

*Monson, Seegmiller , McConnaughey 
& Chen (1989,1990) (Air flow, Re=io s ,io 6 ) 

*Sharma etal (1987) (Axisymmetric 
TAD air flow, Re=io s ) 
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• Earlier Computations: 

*Chen and Sandbom (1986) (K-e and 
curvature-corrected K-e ) 

*Monson, Seegmiller & McConnaughey 

(1989,1990) (mixing-length & K-e with 
curvature-correction) 

*Avva etal (1990) (High Re and Low- 
Re K-e ) 

*Gallardo & Lakshminarayana (1993) 
(curvature-modified K-e ) 

• Agreement in above computations are 
not satisfactory. 
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s/H 


Fig.l Static pressure coefficient on turnaround 
duct inner and outer walls 
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Fig. 2 Skin friction coefficient on turnaround 
duct inner and outer walls 



U/Urn 

Fig.3 (a) Longitudinal velocity in turnaround 
duct, x/H=-4 
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Fig.3(b) Turbulent kinetic energy profile 
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Fig.5(b) Longitudinal velocity in 
turnaround duct, theta=90 deg. 
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Fig.5(c) Turbulent kinetic energy 
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Fig.5(d) Turbulent shear stress 
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Fig.5 (e) Comparison of K and 
0.75*(uu+vv) 
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Fig.7(a) Longitudinal velocity in turnaround 
duct, x/H = 2 (downstream of turn) 
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Fig. 7(b) Turbulent kinetic energy 
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VKI Turbine Nozzle Guide Vane Cascade 


* Measurement by Arts etal (1990) at 
Von Karman Institute 

* M(inlet)=0.15, M(outlet)=0.7 to 
1.11, Re=0.5 - 2 x io 6 , To=415k, T(wall) = 
300k 


* Geometry and grid 



3 . Smoothing: Typical values for 

2nd and 4th order dissipation 
taken as 2% -3% and 3% -4% for the 
turbine cascade computations. 

4 . Table: Computed cases : 


cases 

Mur228 

Mur224 

Mur239 

Pol(bar) 

0.915 

0.909 

3.387 

T 0 1 (K) 

403 

403 

412 

Re, 2 

0.6E + 6 

0.6E+6 

2.1 E+6 

Tu in (%) 

1 

6 

6 


T(wall) = 300k for all these cases. 
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Fig. 9 Blade isentropic Mach number 
distribuition for Mur043 (P ol =1.435 bar, =0.93, 
Tu (inlet) =1%) 
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CONCLUSIONS 


*A two-dimensional Navier-Stokes 
code has been developed using an explicit 
Runge-Kutta method incorporating the 
ARSM model, NLSM model , Chien's LRN 

k-e model and Coakley's LRN q-co model. 

*The surface pressure distributions and 
wake profiles of a transonic turbine cascade 

were predicted well by the k-e , ARSM 

and q-co models. The heat transfer on 
suction surfaces were predicted well by the 

k-e and ARSM models. 

*The heat transfer on pressure side for 
one case (MUR239) was underpredicted. 
This was caused by the underprediction of 
mainstream turbulence level , which 
strongly influences the transition location. 

*The boundary layer code predicts the 
heat transfer on pressure surface well, but it 
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does not capture the transition on suction 
surfaces for all the 3 cases. 

*The wall damping function (c^) in 
Chien's model was modified to yield 
improved prediction for flow under adverse 
pressure gradient. 

*For TAD flow, good predictions have 
been obtained for the surface pressure 
distribution and skin friction coefficients. 
The ARSM model yields better prediction 

than NLSM and k-£ models, for both the 
mean and the turbulence quantities. 

*The near wall "echo" term <j> iiw is not 

correctly modelled for strongly curved 
flows, especially near the concave surface. 

*For more accurate prediction of 
strongly curved TAD flows, it may be 
necessary to use existing RSM, including 
the modeling for the wall region. 
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